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We review basic features of cosmological models with large-scale classical non-Abelian 
Yang-Mills (YM) condensates. There exists a unique SU(2) YM configuration (gen- 
eralizable to larger gauge groups) compatible with homogeneity and isotropy of the 
three-space which is parameterized by a single scalar field. In the past various aspects 
of Einstein- Yang-Mills (EYM) cosmology were discussed in the context of the Early 
Universe. Due to conformal invariance, solvable EYM FRW models exist both on the 
classical and quantum levels. To develop the YM model for dark energy one has to 
find mechanisms of the conformal symmetry breaking. We discuss the Born-Infcld 
generalization and some phenomenological models motivated by quantum corrections 
exploring possibility of transient DE and phantom regimes. 



1 Introduction 

Vector fields as the dark energy candidates were recently proposed by several author ^^ ' 4 ™ 7 ! 
An attractive feature of this suggestion is the fact that vector fields are the basic con- 
stituents of the standard model and its generalizations, while the scalar fields invoked 
in the most popular dark energy models^ remain rather speculative. Formation of non- 
Abelian condensates in cosmology within the context of the standard model and GUT-s 
may have various reasons'^. An inflationary model driven by a vector field was proposed 
by Ford^ long ago. Basically the same ideas lie behind more recent proposals of vector 
dark energy. A major problem with vector fields is that they generically single out spa- 
tial directions leading to anisotropy, unless one is just looking for the source of possible 
anisotropy^l To ensure isotopy one has to introduce suitable multiplets of vector fields. 
A natural framework for this is non-Abelian gauge theory. 

A remarkable fact is that the SU{2) YM field allows for isotropic and homogeneous 
configurations parameterized by a single scalar function. This ansatz is compatible with 
the FRW models of any type, closed, open or spatially flat^ (the compatibility of the 
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closed FRW model with the SU(2) YM, quite obvious due to isomorphism SU(2) ~ 

S 3 , was noticed earlier by several authors D-4 | 15 | 16 |j ^his property remains valid also 

for larger gauge groups containing an embedded SU(2) H ? I 18 | 19 | 20 J21J Remarkably, the 

EYM cosmology is also solvable at the quantum level, the quantum FRW cosmology 

having been discussed in a number of papers ^^^^^^ i An interesting feature the 

EYM quantum cosmology is possibility of tunneling transitions between de Sitter and 

hot FRW cosmologies. An embedding of the EYM cosmology into string theory was also 
discussed [ 2Bl29130| 

Dynamics of classical YM fields in the homogeneous and isotropic FRW cosmology 
is perfectly regular. In the homogeneous but anisotropic models the YM field in con- 
trary often exhibit chaotic 

featuresEMES] 

This behavior is similar to chaotic behavior 
of homogeneous (depending only on time) YM fields in flat space ^3. A very peculiar 
dynamics was observed in the Kantowski-Sachs non-Abelian cosmology, which was inves- 
tigated as an interior solution for non-Abelian black holes^SEZEHEHHD] ^ j g non _ c h ao tic, 
but the approach to singularity is characterized by infinitely growing oscillations. Other 
attempts to apply YM fields in cosmology amount to the idea of growing defects, it was 
applied to Einstein- Yang-Mills sphalerons 1 4 1 1 42 1 m EH] p or more general applications of 
self-gravitating YM fields see Refs.^^1 

Dynamics of the YM fields governed by the standard action quadratic in the field 
strength is conformally invariant (in four dimensions). In the cosmological context this 
means that the equation of state arising in the framework of the FRW cosmology will 
be that of the photon gas p = e/3, therefore one actually deals with a cold substitute 
for the hot gas^l It was even argued that the classical YM condensate could be part 
of the CMB Conformal invariance implies that if the YM condensate was formed 
before inflation, it would be strongly diluted during the de Sitter era. Meanwhile, confor- 
mal invariance might be broken for various reasons, such as string theory modifications, 
quantum corrections or interactions with other fields. Within the context of string theory 
one deals with the Born-Infeld type action for YM fields rather than the Maxwell action, 
which is no more conformally invariant (for discussion and references see®}). The FRW 
cosmology driven by a non-Abelian Born-Infeld action was studied in detail in Ref.^l 
The equation of state in this model interpolates between the string gas equation p = —e/3 
and that of the hot photon gas. The negative pressure may generate the regime of zero de- 
celeration, but it is still insufficient to ensure cosmic acceleration. Further generalizations 
of this model were considered m l 50 | 50 | 51 | 52 | 5 [ y[ OYe general phenomenological models 
with vector fields governed by dynamics of the Born-Infeld type were also discussed in^31. 
Note that in anisotropic cosmology the Born-Infeld Lagrangian typically leads to damping 



off the chaotic oscillationsESEEIHT] 

Here we briefly discuss various aspects of cosmological models with classical large 
scale YM fields which may be helpful to better understand the possibility of the corre- 
sponding dark energy (DE) models. The concrete suggestions in this directions still look 
preliminary, while the idea certainly deserves further investigation. 



2 SU(2) homogeneous and isotropic cosmology 

2. 1 YM configuration 

Space-time symmetry of the Yang-Mills field means that its variation under action 
of the isometry transformations can be compensated by a suitable gauge transformation. 
In this case the field strength F^ u transforms multiplicatively so that the action and the 
stress-tensor remain invariant, implying compatibility with the FRW space-time metric 

ds 2 = a 2 {ri)(drf-dlt). (1) 

Demanding such invariance under spatial translations and SO (3) rotations one finds that 
the unique (up to scaling) SU{2) YM configuration is parameterized by a single scalar 
function of time (we use conformal time) w(rj) and contains both electric and magnetic 
components: 

E? = w6?, B? = (k-w 2 )5?, (2) 

where k = ±1,0 for closed, open and spatially flat cases (for k — 1,0 it was shown 
long ago by Cervero and Jacobs^^l, Henneaux^ and Hosotani^, while derivation for 
all the three cases starting with the Witten ansatz was given by Gal'tsov and Volkov^ 
(for a derivation see also^). The standard YM action quadratic in F" reduces to the 
one-dimensional Lagrangian 

3 (w 2 -(k-w 2 ) 2 ). (3) 
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This looks like the scalar field Lagrangian with the potential. The kinetic term corresponds 
to the (color) electric field contribution, while the potential term is due to the magnetic 
field. During the evolution the electric and magnetic fields exchange the energy, and 
no purely electric or purely magnetic YM fields are possible except the static (unstable) 
configuration w = 0, which is non-vacuum in the closed case k = 1 (the cosmological 
sphaleron). The dependence a~ 4 on the scale factor reflects the conformal nature of the 
YM field. 



2.2 Cold matter for hot Universe 

The standard YM action is conformally invariant, so the above configuration gives rise to 
the equation of state of the photon gas 

p = - with e = - — j \w 2 + (k- w 2 ) 2 } . (4) 

Thus one obtains the hot Universe driven by cold matter 

The YM equations reduce to an equation of motion of a particle in the potential well 

V =(k-w 2 ) 2 J (5) 

which in the closed case (k — 1) is the double well potential. Two minima correspond to 
neighboring topologically distinct vacua. Depending on the w-particle total energy, the 
field classically oscillates either around a single vacuum, or around both of them. In the 
flat and open cases the potential has one minimum at w = with V = for the flat and 
V = 1 for the open cases. The Chern-Simons 3-form 

p 2 / lie \ 

u 3 = — Tr [AAdA-—AAAAA), (6) 



satisfying the equation 
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du 3 = ——TyF A F, (7) 
is non-trivial in the case k = 1. The Chern-Simons number is equal to 

Ncs = J s3 u 3 = j(w + 1) 2 (2 - w). (8) 

The vacuum w = — 1 is topologically trivial: Ncs — 0, the vacuum w = 1 is non-trivial: 
Ncs — 1- Due to the axial anomaly, the evolution of the Chern-Simons density will be 
accompanied by the evolution of the fermion number density for fermions interacting with 
the YM field. 



2.3 Cosmological sphaleron 

In the closed case there exists a particularly simple configuration w = (u>-particle 
sitting at the top of the barrier between two vacua) which in analogy with sphalerons in 
the Weinberg-Salam (WS) theory was called the "cosmological sphaleron"^. I is worth 
noting that the localized particle-like EYM solutions similar to the WS sphalerons exist as 
well^D which are asymptotically flat regular particle-like solutions of the EYM equations 
discovered by Bartnik and McKinnon^Sl. The repulsive YM stresses in these objects are 
compensated by gravity instead of the Higgs field in the WS sphalerons. Creation and 



decay of sphalerons generates a transition of the YM field between topological sectors, and 
it is accompanied by the fermion number non-conservation in presence of fermions^EI] 
The cosmological sphaleron has the topological charge Ncs = 1/2 like the sphaleron in 
the Weinberg-Salam theory. 

The equation of motion of the w-particle 

w = 2w(l- w 2 ) (9) 

is solved indeed by w = 0, this correspond to the total energy 

w 2 + (l-w 2 ) 2 = 1. (10) 

The YM field in this case is purely magnetic. A more general solution with the same 
energy describes rolling down of the w-particle (the sphaleron decay) 

cosh(^/2^) ^ 
Rolling down to the vacuum w — 1 takes an infinite time, while the corresponding full 
cosmological evolution is given by 



a = J—^smr] (12) 

and takes a finite time. Thus the cosmological sphaleron is quasi-stable. This conclusion 
is not modified if a positive cosmological constant is added^. 

2.4 Instantons and wormholes 

An homogeneous and isotropic EYM system has interesting features also in the k = 1 
space of Euclidean signature which is invoked in the path-integral formulation of quan- 
tum gravity. Actually, when \w\ < 1, transitions between two topological sectors can 
be effected via underbarrier tunneling described by instanton and wormhole Euclidean 
solutions. In the Euclidean regime the first integral of the equations of motion reads: 

w 2 - (w 2 - l) 2 = -C, (13) 

where C is an integration constant. Instanton corresponds to C = 0, it describes tunneling 
between the vacua w — ±1. It is a self-dual Euclidean YM configuration for which the 
stress tensor is zero, therefore a conformally flat gravitational field can be added just as 
a background. 

Tunneling solutions at higher excitation levels < C < 1 are not self-dual. In flat 
space-time they are known as a meron (C = 1, the Euclidean counterpart of the cosmo- 
logical sphaleron, N C s = 1/2) and nested merons < C < 1, 1/2 < N C s < 1- The 



energy-momentum tensor of the meron is non-zero, and in flat space this solution is sin- 
gular. When gravity is added, the singularity at the location of a meron expands to a 
wormhole throat, and consequently, the Euclidean topology of the space-tme transforms 
to that of a wormhole. Topological charge of the meron wormholes is zero, the charge of 
the meron being swallowed by the wormhole ^^^^^7 | 23 | 24 | The total action of these 
wormholes diverges because of slow fall-off of the meron field at infinity, so the amplitude 
of creation of baby universes associated with the Euclidean wormholes is zero. However, 
when a positive cosmological constant is added (inflation) the action becomes finite due 
to compactness of the space. Such solutions can be interpreted as describing tunneling 
between the de Sitter space and the hot FRW universe. 

Adding the positive cosmological constant, we will get similar first integrals both for 
the w-particle and the cosmological radius: 

w 2 - (w 2 - l) 2 = -C, a 2 + (Aa 4 /3 - a 2 ) = -C/(e 2 m 2 Pl ). (14) 

Solutions describe independent tunneling of w and a with different periods T w , T a depend- 
ing on the excitation level C. To be wormholes, they must obey a quantization condition 
n w T w = n a T a ^ with two integers. However, for a specific value of the cosmological 
constant 

A = -m 2 Pl , (15) 

it was found^that T a = T w for all C G [0, 1]. In particular, for C = 1 (the meron limit) 
the radius a becomes constant (Euclidean static Einstein Universe). For C ^ 1,0 the 
solutions describe creation of the baby universes (which was invoked in the Coleman's idea 
of the "Big fix"). Remarkably, under the above conditions, the total action (gravitational 
plus YM) is precisely zero^S 

Sym + S gr = 0. (16) 
Thus, the pinching off of baby universes occurs with the unit probability. 

3 Non-Abelian Born-Infeld (NBI) 

Open string theory suggests the following generalization of the Maxwell Lagrangian (ap- 
plicable to any dimensions): 

P 2 



L = ^ (y- det(^ + F^/P) - S^g) , (17) 

P being the critical BI field strength (/? = l/2ira' in string theory). In four dimensions 
this is equivalent to 



P 2 



i+%^- ( ^r )2 - as) 
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In the non-Abelian case the strength tensor F^ u is matrix valued, and the prescription 
(TseytlinS^) is more complicated: the symmetrized trace, which is calculated expanding 
the Lagrangian in powers of F£ v T a (T a being the gauge group generators), then sym- 
metrizing all products of T a involved and only afterwards taking the trace. Symbolically 
this is given by the expression 

T , str = ^Str (j-detig^ + F^/p) - y/=g) , (19) 



but actually this is a useful form if one is able to perform a subsequent resummation. 
Fortunately, this is possible in the closed from for the homogeneous and isotropic SU(2) 
YM field and the metric ds 2 = N 2 dt 2 - a 2 dl\ leading to^H 

l-2£+2^-3^£ w 2 2 (w 2 - k) 2 

A simpler (the ordinary trace) prescription for the NBI Lagrangian consists in summation 
over color indices in the field invariants Ff^ u F atJ-u , F^ u F aflu in the square root form of the 
Lagrangian fllSH . This gives 



■JVaVl - 3K 2 + 3V 2 - 9K 2 V 2 . (21) 



3. 1 NBI cosmology 



Homogeneous and isotropic NBI cosmology with an ordinary trace Lagrangian turns out 
to be completely solvable by separation of variables^. It leads to an interesting equation 
of state: 

P=P^\, (22) 
3(e c + e) 

where e c = /3/47T is the critical energy density, corresponding to vanishing pressure. For 
larger energies the pressure becomes negative, its limiting value being 

p = -e/3. (23) 

This is the equation of state of an ensemble of non-interacting isotropically distributed 
straight Nambu-Goto strings (which indicates on the stringy origin of the NBI Lagrangian). 
In the low-energy limit the YM equation of state p = e/3 is recovered. Thus, the NBI 
FRW cosmology smoothly interpolates between the string gas cosmology and the hot 
Universe. The energy density is 



From the YM (NBI) equation one obtains the following evolution equation for the 
energy density: 

a e + e c 

which can be integrated to give 

a A (e + 2e c )e = const. (26) 

From this relation one can see that the behavior of the NBI field interpolates between 
two patterns: 1) for large energy densities (e ^> e c ) the energy density scales as e ~ a~ 2 ; 
2) for small densities e<e c one has a radiation law e ~ a~ 4 . 

Remarkably, the equation for the scale factor a can be decoupled (g = f3G): 

2ga(a 2 + k) 
a = ~ 2ga? + 3(d 2 + k) 

and admits the first integral 

3 (a 2 + k) 2 + Ago 2 (d 2 + k)=C, (27) 
which allows to draw phase portraits for different k 

• Closed. The only singular point is a = 0, b = 0, (b — a) which is a center with 
the eigenvalues ±iy/6g. Solutions evolve from left to right in the upper half-plane 
as time changes from — oo to oo, and from right to left in the lower half-plane. All 
solutions are of an oscillating type: they start at the singularity (a = 0) and after 
a stage of expansion shrink to another singularity. The global qualitative behavior 
of solutions does not differ substantially from that in the conformally invariant YM 
field model, except near the singularity: 

a(t)=bot-^-t 3 + 0(t% (28) 

where b is a free parameter. Absence of the quadratic term means that the Universe 
starts with zero acceleration in accord with the equation of state p — e/3 at high 
densities. 

• Spatially flat. There is a singular line b = each point of which represents a 
solution for an empty space (Minkowski spacetime). This set is degenerate, and 
there are no solutions that reach this curve for finite values of a. All solutions in 
the upper half-plane after initial singularity expand infinitely. A remarkable fact is 
that for this case one can write an exact solution for a in an implicit form: 

A^g (t -t ) = Vs(n- arctan IT 1 + tt/2) , (29) 
where VI = \[2aj \J Va 4 + C — a 2 . The metric singularity is reached at t — t . 



• Open. Physical domain is b < — 1, b > 1. There is a center at a = 0, 6 = with the 
eigenvalues ±iy/6g, but it lies outside the boundary of the physical region. Other 
singular points are (a = 0, b = ±1). These points are degenerate and cannot be 
reached from any point lying in the physically allowed domain of the phase plane. 
The only solutions which start from them are the separatrices b = ±1 that represent 
(part of) the flat Minkowski spacetime in special coordinates. One can easily see 
that all solutions in the upper part of the physical domain a > 1 start from the 
singularity and then move to a — > oo, a — > 1. 

3.2 NBI universe with A 

Cosmological constant gives rise to new types of solutions including de Sitter-like, which 
are non-singular. However, in the case when there is a singularity, its structure is de- 
termined by the leading term in equation of state, namely, by the BI pressure, and is 
unaffected by A. The generic solution near the singularity starts as 

a(t) = b t - ^j^b t 3 + 0(t 5 ), (30) 

with zero acceleration. The equation for the scale factor can be separated again and 
admits the first integral 

3 (V -±a 2 + k} + Aga 2 (b 2 - K 2 + kj=C, (31) 

which gives phase portraits for different k. The allowed domain now is 

9 A , 
a - —a + k >0 

For k — 1 and a negative cosmological constant the only singular point is (a = 0, b = 0), 
which is a center with eigenvalues £ = ±3i^3(2g — A). All phase space trajectories are 
deformed circles, and all solutions start from the singularity and reach the singularity in 
the future. For A > the allowed domain is bounded by a hyperbola b 2 — |-a 2 > — 1 
which represents the de Sitter space. Behavior of the solutions depends on whether A is 
smaller or greater than the quantity 2g. 

• Closed with A. Physical domains are bounded by hyperbolas. The point a = 0, b = 
remains a center. Another pair of singular (saddle) points a = ±y ^^~^ -, 6 = 
lie within the allowed domain and correspond to the static Einstein universe. Enter- 
ing them separatrices correspond to either a solution developing from the singularity 
into the static universe, or one rolling down from an infinite radius to the static uni- 
verse. They divide the phase space into the domains containing different types of 



generic solutions. Near the origin, all solutions evolve from an initial to a final sin- 
gularity. In the region near the physical boundary, the solutions are non-singular 
and evolve for an infinite time first shrinking to some finite value of the scale factor 
and then ever expanding. Generic solutions of the third type possess a singularity 
but are non-periodic: the universe expands forever. 

• A > 2g With increasing A, the saddle point approaches the origin and finally, when 
A > 2g, swallows it. The character of the singular point (a = 0, b = 0) changes — it 
becomes a saddle point. The separatrices entering it are solutions of an inflationary 
type. They divide all generic solutions into two classes: nonsingular de Sitter-like 
and ever-expanding solutions with an initial singularity. 

• Open and flat with A. For k = the only singular point is the origin (a = 0, b = 
0) which now is degenerate. For A < the allowed domain is the whole plane. All 
solutions are of an oscillatory type evolving from the initial to the final singularity 
for a finite time. 

For a positive cosmological constant the allowed domain consists of the upper and 
lower parts of a cone b 2 — ^a 2 > whose boundary corresponds to the (part of) de 
Sitter space described in the inflating coordinates. All solutions that lie within the 
allowed domain have an initial singularity and are ever expanding in the future. 

For k = — 1 the allowed domain lies outside of either an ellipse (A < 0), or a 
hyperbola (A > 0). In this case two types of singular points — the origin and a = 
± y A(4g~j f^? b = (when they are real) lie outside the allowed domain. Another 
pair of singular points is (a = 0,b = ±1). These are degenerate (as in case with 
zero cosmological constant). Since they lie on the boundary of the allowed domain, 
the only physical solutions approaching them correspond to zero energy of the BI 
field ((anti)-deSitter). 

All generic solutions within the allowed domain possess a singularity and are either 
oscillating (for A < 0) or ever expanding (for A > 0). 

3.3 Dynamics of YM field 

The NBI equation for w can be integrated once, giving a first order equation 

(32) 



• 2 4 



a\-{k-w 2 ) 2 a 4 + 3c4' 

whose solution is oscillatory and is given by Jacobi elliptic functions. The main difference 
with the ordinary EYM cosmology relates to small values of a. One can see that near 



the singularity (a — > 0) the YM oscillations in the NBI case slow down, while in the 
ordinary YM cosmology the frequency remains constant in the conformal gauge. Near 
the singularity 

w = w + -^ W r1^ + O{1*), (33) 
6(fc + 6g) 



where a = ±y3(k + 6q) 2 — Ag 2 (k — Wq) 2 , w is a free parameter, b is a metric expansion 
parameter. 

3.4 NBI on the brane 

Replacement of the standard YM Lagrangian by the Born-Infeld one breaks conformal 
symmetry, providing deviation from the hot equation of state and creating negative pres- 
sure. Surprisingly enough, putting the same NBI theory into the RS2 framework gives rise 
to an exact restoration of the conformal symmetry by the brane non-linear corrections^^. 
Choosing the ordinary trace action 

S = XTrj yJ-detig^ + F^/P) d 4 x - k 2 J (R 5 + 2A,) v ^f 5 d 5 x, (34) 

where the brane tension A plays a role of the BI critical energy density, one obtains the 
constraint equation 

(;M A+ s< A+E > 2+ 7-?- < 35 > 

where S is the integration constant corresponding to the bulk Weyl tensor projection 
( "dark radiation" ) and, as usual, 

A 4 = ^ 2 (A + ^ 2 A 2 ), G(4) = ^. (36) 

The energy density in this model scales as 

e = A (^(1 + C/a* - l) , (37) 

where C is the integration constant. Surprisingly, the constraint equation comes back to 
that of the YM conformally symmetric cosmology 



a\ 2 8nG(d\ . C 

cr a 



>A) ^+~--^ (38) 



where the constant C = £ + k 4 A 2 C/36 includes contributions from both the "dark radia- 
tion" and the YM energy density. 



3. 5 Smoothing YM chaos with NBI 

In flat space the homogeneous YM dynamics is typically chaotic^. When gravity is 
switched on one deals with cosmological solutions. The occurrence of chaos is then related 
to additional symmetries: in the isotropic cosmology case, as we have seen, the evolution is 
perfectly regular, but once an isotropy is broken, the chaotic behavior is typically manifest 
^ see |32l33] £ or Bjanchj t and^for other anisotropic Bianchi types). When the quadratic 
YM action is changed to the NBI form, it turns out that YM evolution becomes more 
regular. This may be interpreted as smoothing stringy effect on classical chaos 
For the Bianchi I metric 



with N, b and c depending on time, the corresponding YM field is parameterized two 
functions u, v of time: A = T\udx + T 2 udy + T 3 vdz. The Poincare sections signal the 
chaos-order transition at /3 cr = 0.317 (regular from the strong NBI side (3 < /3 CT ). 

4 Conformal symmetry breaking and DE 

Conformal symmetry breaking in NBI theory demonstrates the occurrence of the negative 
pressure, but its extremal value p = — e/3 is still insufficient for DE. Meanwhile, a stronger 
violation of conformal symmetry may provide an equation of state with e ~ —1. Such 
violation can be of different nature: 

• Quantum corrections, 

• Non-minimal coupling to gravity, 

• Dilaton and other coupled scalar fields including Higgs, 

• String theory corrections. 

Here we just explore some model Lagrangians to see the necessary conditions for DE. 
Assuming the Lagrangian to be an arbitrary function L{T ', Q) of invariants 



ds 2 = N 2 dt 2 - b 2 (dx 2 + dy 2 ) - c 2 dz 2 , 



(39) 



T 



F"F a ^/2 and Q 




(40) 



one finds for the pressure and the energy density (conformal time): 




(41) 



For a simple estimate consider the power-low dependence: 



L ~ T{Tl^y-\ (42) 

where /i has the dimension of mass. Then in the electric (kinetic) dominance regime one 
obtains 

W = V -= 3 ~ 2U (43) 
e 3(2i/ -1) K } 

For certain v this quantity may be arbitrarily close to W = — 1 or even less. An electric 

phantom regime is thus possible. 

In the magnetic (potential) dominance regime one obtains 

W = 4i//3 - 1. (44) 

The value W = — 1 can not be reached, but an admissible DE regime is also possible. 
These regimes are transient since during the evolution the electric part transforms to the 
magnetic and vice versa. 

Another plausible form of the lagrangian (suggested by quantum corrections) is loga- 
rithmic^: 

L~^ln(^//i 2 ). (45) 

Then the energy density is 

e = 3 (T[\n(F/iJ 2 ) + 2] + V ln(^//i 2 )) a" 4 , (46) 

and the equation of state is 

W = P -= — + V) ln ^ 2 ) + 2{2V ~ T ) (A 7 ) 
e 3(T + V) MF/fi 2 ) + 6T ' 1 1 

where T = w 2 , V = (k — w 2 ) 2 . It is easy to see that W ~ — 1 for ln(jF/ fj 2 ) ~ —1. In this 
case the DE regime is transient. The phantom regime is also possible. 

Thus, the DE conditions for the homogeneous and isotropic YM field can arise indeed 
as a result of sufficiently strong breaking of the conformal symmetry. 



5 Outlook 

Attractive features of the hypothesis of the non-Abelian condensates as dark energy are 
clear physical meaning of the Yang-Mills field, the existence of the natural isotropic and 
homogeneous configuration removing the problem of isotropy of Abelian vector field cos- 
mological models, as well as some non-trivial effects related to the YM topology. The 
major yet unsolved problem is the mechanism of the conformal symmetry breaking which 



is necessary to ensure an accelerating expansion. Various ideas for this were suggested, but 
further work is needed to support or reject the hypothesis. Also, physical mechanisms of 
generation of non-Abelian condensates in cosmology require more detailed investigation. 
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